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Abstract 

We derive high-order threshold corrections for top quark production in hadronic col- 
lisions from resummation calculations. We present analytical expressions for the cross 
section through next-to-next-to-next-to-next-to-leading order (N 4 LO) and next-to-next- 
to-leading logarithmic accuracy. Special attention is paid to the role of subleading loga- 
rithms and how they relate to the convergence of the perturbation series and differences 
between various resummation prescriptions. It is shown that care must be taken to avoid 
unphysical terms in the expansions. Numerical results are presented for top quark produc- 
tion at the Tevatron. We find sizeable increases to the total cross section and differential 
distributions and a dramatic reduction of the factorization scale dependence relative to 
next-to-leading order. 
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1 Introduction 

The top quark production cross section at the Tevatron receives significant contributions from 
the threshold region, where there is limited phase space for the emission of real gluons. The 
incomplete cancellation of infrared divergences between real and virtual graphs produces finite, 
but large, logarithmic corrections in the form of "plus" distributions. These corrections can be 
resummed to all orders in the perturbative expansion. 

The need for threshold resummations was recognized over a decade ago for the Drell-Yan 
cross section [|. The resummation of the leading logarithms for heavy quark production 
0. [H H relies heavily on this work since these logarithms are universal between electroweak and 
QCD induced cross sections. To go beyond leading logarithms one has to take into account the 
complex color structures of QCD cross sections || [7|. Resummation for heavy quark production 
at next-to-leading logarithmic (NLL) accuracy was first presented in Refs. [0, [7| and then in 
Refs. i, 0, [Hi 0, |TJ. For a review see Refs. @|. 



Threshold resummation follows from the factorization properties of the cross section [|T5|]. 
separating long- and short-distance physics. The hadronic cross section is written as a convo- 
lution of non-perturbative parton distribution functions with the perturbative partonic cross 
section. This convolution becomes a simple product if one takes moments of the cross section. 
The perturbative cross section still has sensitivity to soft-gluon dynamics and is a smooth func- 
tion only away from the edges of partonic phase space. In general it includes "plus" distributions 
with respect to a variable that measures distance from partonic threshold. It is these singular 
distributions that can be resummed to all orders in perturbation theory. This is achieved by 
first refactorizing the cross section into new functions which absorb the universal collinear 
singularities associated with the incoming partons, and a function that describes non-collinear 
soft gluon emission. Resummation is explicitly derived in moment space from the renormal- 
ization group properties of these functions [0, [TT], [16], |17| . To obtain the physical resummed 



cross section the moment space results must be inverted back to momentum space. There have 
been a few proposals or "prescriptions" for the best way to do the inversion or otherwise use 
the resummation [[3], [|, [5]. Numerically, the choice of prescription can have a significant effect. 
Alternatively, the resummed cross section can be used as a generator of perturbation theory. 
At fixed order no prescription is necessary to invert the moment-space results. This is the 
approach taken in this paper. 

The aim of this paper is two-fold. First, to derive the expansion of the resummed cross sec- 
tion through next-to-next-to-next-to-next-to-leading order (N 4 LO) and next-to-next-to-leading 
logarithmic (NNLL) accuracy; and second, to investigate and assess the importance of sublead- 
ing logarithms (beyond NNLL). We calculate the numerical effect of these higher orders and 
subleading logarithms for top quark production at the Tevatron. As we will see, this discus- 
sion intimately relates to the convergence of the perturbation series as well as to the choice 
of a resummation prescription, over which there have been different viewpoints || |18| . We 
will show that one has to be careful to avoid unphysical terms in the fixed-order expansions. 
These considerations apply to QCD hard scattering cross sections in general. We find large 
corrections for top quark production at the Tevatron and a dramatically reduced factorization 
scale dependence. Our formalism allows us to make predictions for both total and differential 
cross sections. 



In a companion paper [19| the next-to-leading order (NLO) and next-to-next-to-leading or- 



der (NNLO) threshold corrections to heavy quark production at NNLL accuracy are also studied 
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and the difference between single-particle inclusive (1PI) and pair inclusive (PIM) kinematics 
is highlighted. Similar NNLO expansions have recently been presented for electroweak-boson 



20], [211 , direct photon |22|, [23| , and jet [gj hadroproduction. 

In the next section we briefly review the resummation formalism and give the expression 
for the resummed cross section in moment space. In Section 3 we present NLO and NNLO 
expansions at NNLL accuracy in single-particle inclusive kinematics and we discuss in detail the 
role of subleading logarithms. We also present numerical results for the top total cross section 
and transverse momentum distributions at the Tevatron. In Section 4 we present the next-to- 
next-to-next-to-leading order (N 3 LO) corrections at NNLL accuracy, and again we study the 
role of subleading logarithms at that order. In Section 5 we discuss N 4 LO and higher-order 
corrections. In Appendix A we list several formulas for the Mellin transforms that are used 
in the calculations. In Appendix B we collect some formulas for the integrations involved in 
the calculation of the hadronic total and differential cross sections. In Appendix C we present 
results for the NLO and NNLO expansions in PIM kinematics. 



2 Resummed cross section 

We begin by briefly reviewing the resummation formalism of Refs. || [?[ [TTJ. We will mostly 
discuss the single-particle inclusive cross section but will make a few comments on pair inclusive 
kinematics where appropriate. The differences between the two kinematics in the expression 
for the resummed cross section are minimal. 

The factorized single-heavy quark inclusive cross section for hadron-hadron collisions, 

h A (pA) + h B ( PB ) -> Qfa) + X[Q] , (2.1) 

where the Tt's are the colliding hadrons, Q is the produced heavy quark, and X represents the 
additional partons in the final state including the heavy antiquark, takes the form of a convo- 
lution of the perturbative short- distance cross section a with the universal parton distribution 
functions 6: 



&h A h B ^QX 



dx a dx b <pf a /h A (x a ,n 2 F ) (f>f b /h B (xb, /4) Vfaf b ->Qx (s4,t 1 ,u 1 ,m 2 ,/i 2 F ,a s (f/ R )^ , 

(2.2) 

where fip and llr are the factorization and renormalization scales, respectively. Note that here 
and in the following o and a can denote either the total cross section or any relevant differential 
cross section. The parton processes involved at the Born level are 

?(Po) + q(pb) -> Q(Pi) + Q(P2) , 

g(p a ) + g(p b ) - Q(pi) + Q(p2). (2.3) 

The qq channel is dominant for top production at the Tevatron and contributes over 90% of 
the cross section at the Born level. The partonic invariants in Eq. (2.2) are defined by 



(Pa + Pb) , t 1 = (p a -pi) -m, ux = (p b -pi) - m , (2.4) 



with m the heavy quark mass, while S4 = s + £1 + U\ measures the distance from threshold; at 
threshold S4 = 0. 



4 



By taking moments and replacing the incoming hadrons by partons, the convolution in Eq. 
2]) simplifies to a product 0, [H] , 

*Uh-+Qx(N) = h a /fa(N a ) fa /fb (N b ) *f afb ^Qx(N) . (2.5) 

The moments are defined by &(N) = J(dsi/s) e~ Ns4 ^ s a(s4), with A" the moment variable, and 
similarly for the 0j's with respect to Xi. The definition of iVj depends on the kinematics. For 
1PI kinematics, N a = N(-u/s) and N b = N(-t/s) WU. For PIM kinematics, N a = N b = N 



The short- distance perturbative cross section a still displays sensitivity to soft gluon emis- 
sion. The incomplete cancellation between graphs with gluon emission and virtual gluon correc- 
tions manifests itself in the occurence in a of "plus" distributions which are singular at = 0, 
the partonic threshold. At nth order in a s (beyond the Born term), these distributions are of 
the form [(ln fc (s4/m 2 ))/ £4]+, k < 2n — 1. The leading logarithms correspond to k = In — 1, NLL 
to k = 2ra-2, NNLL to k = 2n-3 and so on. Under moments [(\n 2n ~ 1 (s 4 /m 2 ))/s 4 \ + -> \n 2n N 
(see Appendix A), and our goal becomes to resum logarithms of N. 

To separate these soft gluon effects from the hard scattering, a refactorization is introduced 

&Uh-+QxW = $f«/fM $ Mh (N h ) H{f S f jf(m/(N^ F )) , (2.6) 

where the r/>'s are center-of-mass parton distributions |j] that absorb the universal collinear 
singularities associated with the initial-state partons, and S is the soft-gluon function that 
describes non-collinear soft gluon emission. The mass of the heavy quarks protects the final 
state from collinear singularities. H is the hard-scattering function and is free of soft-gluon 
effects and thus independent of N. H and S are matrices in the space of color exchanges (I, J 
are color indices) and differ for each partonic process. 

Using Eqs. fl2.5p and ( |2.6p to solve for the perturbative cross section a, we then have 
a(N) = (i'f a /f a ipf b /fj(^f a /f a (t)f b /f b )) Tt[HS}. After resumming the N-dependence in ijj/(f) and 
S |], 0, |TT|| , we obtain the resummed heavy quark cross section at NLL accuracy^ in moment 
space: 

dp' 



°fah^Qx(N) = exp 



exp 



x exp 



x Pexp 



V* dp! 



E ifa) {N a ) + E ( - fb) (N b ) 

Tr{H fafb (a s (^L 
S fafb (l,a s (m 2 /N 2 



4( 7a (a s ( yU ' 2 ))+ 7b (a s (/i' 2 )) 



II %(3{a 8 ^ 2 )) 



Pexp 



Jm fi 



(2.7) 



This expression is actually valid for both 1PI and PIM kinematics with appropriate definitions 
for N a and A^ as discussed previously. 

The first exponent in Eq. (|2.7| ) resums the A^-dependence of the ratio and is 

given in the MS scheme by 

^-1 _ l f r & dfi' 2 



dz- 



1 - z 




(l-z) 2 s fj, 



12 



1 

+ 2 K 



(fi) 



(2.i 



■"■Note that although we formally have NLL accuracy in the resummed exponent, after matching with the 
exact NLO cross section we can determine all NNLL terms in the finite-order expansions. 
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At next-to-leading order accuracy in IniV, we need A^> at two-loops, A^\a s ) = Cf[a s /ir 
+ (a s /n) 2 K/2], and = 2C f (a s /n)(l - ln(2i//)). Here C f = C F = (N 2 - 1)/(2JV C ) for an 
incoming quark, and Cf — Ca = N c for an incoming gluon, with N c the number of colors, while 
K = Ca (67/18 — 7r 2 /6) — 5n//9, where rif is the number of quark flavors. The Ui terms are 

gauge dependent. They are defined by i>i = ($ ■ n) 2 /\n\ 2 , where (3i = pi^2/s are the particle 
velocities and n is the axial gauge vector. We note that all gauge dependence cancels out in 
the cross section. 

In the DIS scheme, which is usually only applied to quarks, the first exponent in Eq. ((27 
is given in terms of the MS result as 



E^\N l )\ ms = E^\N t )\ m - \~dz *\ / { f^A® [a s (Xs)]+B^ [a s ((l - z)s)} 




(2-9) 



where B^(a s ) = -(3C F /4)(a e /ir). 

The 7 a are anomalous dimensions of the fields if) and are given at one loop by 7 q = 
(a,j/7r)(3CV/4) and 7 9 = (a s /7r)(/3o/4) for quarks and gluons, respectively. The (3 function 
is given by (3{a s ) = \id\ngjd\i = — /3oa s /(47r) + with flo = (IICa — 2n/)/3. 

The trace appearing in the resummed expression is taken in color space. The symbols P 
and P denote path ordering in the same sense as the variable \J and against it, respectively. 
The evolution of the soft function from scale m/N to m follows from its renormalization group 
properties and is given in terms of the soft anomalous dimension matrix Tg 0. For the 
determination of T $ an appropriate choice of color basis has to be made. For the qq channel 
we use an s-channel singlet-octet basis, while for the gg channel we use a basis consisting of 
three color tensors . The soft anomalous dimension matrix, evaluated through the calculation 
of one-loop eikonal vertex corrections, has been presented for the partonic processes in heavy 
quark production in Refs. || 0, 0. In the color bases that we use, the soft matrices, S, 
are diagonal for both partonic channels, and the hard matrix for the qq channel has only one 
non-zero element. At lowest order, the trace of the product of the hard and soft matrices 
reproduces the Born cross section in each partonic channel. We also note that the matrices 
are not diagonal in the color bases that we use. If we perform a diagonalization so that 
the matrices do become diagonal, then the path-ordered exponentials of matrices in the 
resummed expression reduce to simple exponentials; however, this diagonalization procedure is 
complicated in practice |T3fl . 



The integrations over z in the exponents of Eqs. (277), ( |2.8| ) run over the region where 
the running coupling constant a s diverges. The prescriptions of Refs. 0, £|, ||] have been 
proposed to avoid these soft gluon divergences in the resummed cross section. However, if we 
expand the exponents in the resummed cross section at fixed order in a s and invert back to 
momentum space using the equations in Appendix A, no divergences are encountered and thus 
no prescription is required. In addition to avoiding the necessity for a resummation prescription, 
a finite-order expansion bypasses the need for the diagonalization procedure that we mentioned 
above, as well. 



3 NLO and NNLO threshold corrections 



In this section we expand the resummed cross section to next-to- leading and next-to-next-to- 
leading orders. In the following stands for the nth-order differential corrections. Nominally, 
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it denotes s 2 d 2 a^ /(dtidui) but it can also denote any other relevant differential cross section, 
such as cPcr*™ - ) / (dp^dsi), with p F the transverse momentum, or d 2 a^ /(dyds^), with y the 
rapidity, or d 2 a^ / (dpxdy), with appropriate Jacobians inserted into the definition of the Born 
term, a B , and the function -Bqed for the gg channel in the expressions below. 



3.1 NLO threshold corrections 



We first expand the resummed cross section to next-to-leading order in 1PI kinematics. These 
expansions are already known for both 1PI and PIM kinematics |7], [TT], |T3| |. 

For the qq channel in the MS scheme, the full next-to-leading-order threshold corrections 

are 



a 



MS(1)_ 
qq^QQ 



's 4 ,m 2 ,S,ti,Ui,/J, F ,fJ, R ) = CT 



+ 



1 

Si 



qq^QQ 



2ReT' q 2 q 2 -2C F + 2C F hi 




. x( \ (l)qqS+V 



(3.1) 



where ^jjg S+V denotes the soft plus virtual (S + V) 8(s±) terms in the NLO cross section that 
can be obtained from Eq. (4.7) in Ref. |25|] (with t\ and U\ interchanged because of different 
definitions in that reference). Also 



Rer™ 



C f 



41n(^ 



ReL, 



Ca 
2 



-31n('^ 



In (^] + ReL 



(3.2) 



is obtained from the real part of the one-loop soft anomalous dimension matrix element T% 
after dropping all gauge dependent terms and an overall coefficient a s /n. Here Lg = (1 — 

2m 2 /s)/j3 ■ [ln((l — (3) /{I + (3)) + m], with (3 = Jl — Am? /s, is the velocity-dependent eikonal 
function. The Born term is given by 



Oqq-^QQ = ^ a]^) K qq N C C F 



tj + uj 



2m 2 



(3.3) 



where K q q = N c 2 is a color average factor. In our calculation for the expansion, we used the 
result for the soft matrix at lowest order in Eq. (A5) of Ref. The lowest-order hard matrix 

has only one non-zero element, given by ifff - ^^ = [2 '/ '(^ T cC 'f)]o"^^qq- 

To be sure, the expansion of the NLL resummed cross section does not give all the S(s4) 
terms, only those 5(s4) terms involving the scale; these terms are 0"^^Qg(a s (/i|.)/7r)(5(s4)[(— 3/2 
+ \n(t 1 Ui/m 4 ))C F \ia(fi F /m 2 ) + (fio/2) ln(/^/m 2 )]. The rest are obtained by simply matching 
with the NLO cross section in p5[ . Thus we obtain all the S + V terms at NLO. As shown 
120, 



in Refs. 



these terms dominate the cross section and are an excellent approximation at 



the partonic level to the exact NLO cross section close to threshold and even quite far from it. 
We note that the exact NLO cross section is the sum of the full S + V terms and hard gluon 
corrections; the latter are not taken into account by threshold studies and vanish at threshold. 
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In the DIS scheme, the corresponding result for the qq channel is 



a 



DIS (1)_ 
qq^QQ 



(s 4 ,m 2 ,s,t l ,u 1 ,fi F ,fi R ) 



a 



« s (/4) 



qq^QQ' 



2C f 



TV 



ln(s4/m 2 



+ 



1 

S4 



2 



2ReT' g 2 l - + C F In 



2C F \n 



S4 

V m 2 



where cr^g S+V can be obtained from Eq. (4.14) in |25 . 

For the gg channel in the MS scheme the NLO threshold corrections are given by 



-MS(1) / 2 + 

a gg^QQ\ S ^ m , S , t X , U ly fl F , fi R 

+ ^s^ R ) K 99 B QED 



a 



B 

gg^QQ' 



a s (jJL 



AC, 



Nc{N ( 



1 

2 Utt, 



ln(s4/m 2 ) 



s 4 



2C A ln 



m 2 



-C F + ^ ) R<>Z , 



(N 2 - 1) 
+ 1 \ T ) (C F -C A )ReL p 



— + ^f(iV 2 



AT, 



(3.4) 



1 

S4 



(3.5) 



where = (A" 2 — 1) 2 is a color average factor, and 



m 2 s 



_ t X Mi 4^5 / 

#qed h — + 1 - 



(3.6) 



Here, &^§ 9S+V again denotes the soft plus virtual <5(s 4 ) terms in the NLO cross section. These 
terms are given by Eq. (6.19) in Ref. [p8[| ; we note, however, that in that reference the scale was 
set equal to m, therefore in addition to those terms we have to include in ct^ ssS+V the terms 



a 



B 



gg ^ QQK a s (/j, 2 R )/n)[C A \n(t 1 Ui/m' 1 ) ln(/4/m 2 ) + (f3 /2) \n(fi 2 R / fx 2 F )}. The Born term is given by 



a 9 9 ^QQ = ^<{^R)K gg N c C F 



C F -C 



A — Y 

s 2 



B, 



QED 



(3.7) 



In our calculation for the expansion we used the result for the soft matrix at lowest order in 
Eq. (C3) of Ref. [p4 . The lowest-order hard matrix has the form of Eq. (C6) of Ref. [24 



with independent elements E{V QQ = 7va 2 s B QED K gg /(2N 2 ), E{V QQ = N C H??~* QQ (if - uj) / s 2 , 

and H 9 3 r QQ = N 2 H(r QQ (l - 4:t lUl /s 2 ). 

Again as shown in Refs. [T2^, E7| these corrections dominate the cross section near threshold 
and are an excellent approximation at the partonic level to the exact NLO cross section. For 
NLO expansions in PIM kinematics see Appendix C. 



3.2 NNLO-NNLL threshold corrections for qq -> QQ 

Next we derive the NNLO-NNLL threshold corrections from the two-loop expansion of the 
resummed cross section. 



s 



For the qq channel in the MS scheme these corrections 



are 



.MS (2) / 2 j \ B 



qq^QQ 



2\\ A 



IT 



8C* 



ln 3 (s 4 /m 2 ) 
s 4 



+ 



+ 



ln 2 (s 4 /m 2 ) 
s 4 

ln(s 4 /m 2 ) 
s 4 







- + 


4 



-PoC F + 12C F f ReT'g -C F + C F In (^^^ 
Rer'g - C F - C F In f ^i) - C F In (4) 



-C F \n ^ 



m z 



+ 4r'?lr'g - /% 



Rel^g - CV - C F In 



Vsm 2 / 



m z 



+ 2C F K - 16C 2 C| + AC F cgf s " 
where d^ 5S+V is defined by 



-V 



} + o 



1 

Ls 4 



(l)qqS+V _ «£ B _(1)WS+V 
°MS ~~ ^ °<?9^QQ L MS 



7T " ■ ■ 

and the off-diagonal elements of the soft anomalous dimension 

a s /Tr) are 



(3.8) 
(3.9) 

matrix (dropping an overall 



r?! = 21n(^ 



V ti J ' " Ca V h J ^ 
We are able to derive all the NNLL terms by matching with the S + V terms in the NLO cross 
section, Eq. (|3.9|). 

In addition, we can derive at NNLL accuracy the following [l/s 4 ]+ and <5(s 4 ) terms involving 
logarithms of the factorization and renormalization scales: 

1 



°qq^QQ 1 



7T 



S4 



3 
-( 

2 



+ 

!0,Aln(4)ln(4)+ln(4 
! \m A J \m A J \m A 

+ C ,(2 h (^)-3)fRaT 



(7p f 3 - 2 In 



-2C F j g> * - C F A' + 8C F < 2 



'22 — CV ~l~ CV In 



Rer'g - C F - C F In 



m 2 s 



1 B 
"I" Vqq^QQ 



TV 



{in 2 f 4 



2C 2 C2-^lnf^+4^ 



+ ^0 In f4) ^ (4 



16' 



In 



m 4 



3 

~~ 2 



^C 2 -^| In 



m 4 



(3.11) 



where T^ qq is obtained from c^" S+V in Eq. (|3.9| ) by dropping all scale terms \n(/j lFl 

ln(// fl /m) in c^"" S+V . We note that at NNLL accuracy we derive all [l/s 4 ] + scale tei 
in the 5(s 4 ) coefficient we can only determine quadratic terms in the scale logarithms. 
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qq->tt S l/2 =1.8TeV 

12 k 



, I I 1 1 | ! 
\ 


i i i 1 i i i i 1 i i i i 1 i i i i 

NNLO-NNLL upper curves 




NLO lower curves 


x xx 

— \ ^"'x 

v Xx 

x \- 


m/2 < \i < 2m 







150 160 170 180 190 200 

Top mass (GeV/c 2 ) 

Fig. 1. Top quark production at the Tevatron with y/S = 1.8 TeV for the qq — > ti channel in 
the MS scheme. Plotted are the exact NLO cross section for fi = m (lower solid line), m/2 and 
2m (upper and lower dashed lines), and the NNLL-NNLO cross section for n = m (upper solid 
line), m/2 and 2m (upper and lower dotted lines). 
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10 



qq->tt m=175GeV/c" S l/z =1.8TeV 



NNLO-NNLL 
NLO exact 




2 



(a./ m 

Fig. 2. The scale dependence of the cross section for qq 
with VS = 1.8 TeV and m = 175 GeV/c 2 . 
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tt in the MS scheme at the Tevatron 



Since the qq channel is dominant for top quark production at the Tevatron, it is worthwhile 
to show some numerical results for the hadronic cross section in that channel alone before 
presenting the full cross section (for details of the hadronic calculation see Appendix B). In 
Fig. 1 we plot the NNLO-NNLL top quark cross section qq — > ti at the Tevatron with = 1.8 



TeV together with the exact NLO cross section [g5|, |29| as a function of the top mass; we have 
set fi = [Ap = fiR. Here and in the rest of the paper we use the CTEQ5M parton densities PU 
when calculating MS results. We note the significant increase of the cross section at NNLO 
along with the dramatic reduction in scale variation between m/2 and 2m. This reduction is 
also evident for a wide range of scale choices in Fig. 2. The NNLO cross section is larger than 
at NLO and relatively flat with respect to scale variations. 

In Table 1 we present detailed numerical results for the qq corrections to the top quark 
cross section at the Tevatron with m = 175 GeV/c 2 through NNLO at both NLL and NNLL 
accuracy. Some of these numbers will be useful in our discussion in the next subsection. 

The NNLO-NNLL threshold corrections in the DIS scheme are 



a 



DIS (2)_ 
qq^QQ 



a 



B 

qq^QQ 



+ 



+ 



ln 2 (s 4 /m 2 ) 
s 4 

ln(s 4 /m 2 ) 
s 4 




C F + 6C F ReT'H - + In 



C_ F 
4 







- + 


4 



ReT'll - 



C f 



— In 



C F ln 
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qq — > ti 


fi — m 


H = m/2 


H = 2m 


Born 


3.81 


5.30 


2.87 


NLO-exact 


0.81 


-0.53 


1.37 


NLO-NLL 


1.31 


0.03 


1.81 


NLO-full S + V 


1.26 


-0.06 


1.78 


NNLO-NLL 


1.01 


0.67 


1.29 


NNLO-NNLL 


0.80 


0.62 


1.22 



Table 1: The MS corrections for top quark production in the qq channel in pb for pp collisions 
with yS = 1.8 TeV and m = 175 GeV/c 2 . Here ji — jJt, F — [Mr. 



+ 4ir 2 r™ - a 



Rer' 22 -^ F -^ln(^ 

4 V mr 

1 

s 4 



-C F ln 



+ CpK - 4( 2 C 2 F + 2C F c b 1 I ) s ?9 " S+V ] } + O 



(3.12) 



where Cuj S 99S+V is defined in analogy to Eq. ( |3.9| ). As for the MS corrections, we can also derive 
additional [1/54]+ and 5(s±) terms involving the scale in the DIS scheme. These terms are 



a 



B 

qq^QQ 



71 



In 



2 I Vf_ 

m 2 



/4 



/4 



^CWfcln In ^ +ln ^ 



2 

+ C F 



C F 3 - 2 In 



2C F rS ) s <?9 "-C F K + 4CX2 

2 \ ' 



+ 



In 



t x u x 



rrr 



2Rer ,w " 



22 



— + CV In 

2 W lUl 



Rer'l-^-^ln 
22 4 2 



B 

qq^QQ 



' a s (jJL 



71 



00 



-2C F ( 2 -^-C F ln 



l4 



5(s 4 

m 4 
'fit 



In 



2 / Afp 

m 2 
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where T$ s 9 ^ is defined in analogy to its MS counterpart. 

At the Tevatron, with y/S = 1.8 TeV and ji F — jiR — m — 175 GeV/c 2 , and with the 
CTEQ5D parton densities |3(|, the exact NLO cross section in the DIS scheme is 4.60 pb and 
the NNLO-NNLL corrections provide an additional 0.30 pb. The NNLO corrections in the 
DIS scheme are much smaller than the corresponding MS corrections; that follows from the 
definition of the two schemes. 

Results for the NNLO expansion of the resummed cross section in PIM kinematics in both 
the MS and DIS schemes are presented in Appendix C (see also Ref. f!3]|). The differences be- 
tween the expansions in the two different kinematics are in the extra terms involving ln(tiMi/m 4 ) 
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for the 1PI fixed-order expansions relative to the PIM expansions and in the matching terms 
needed to reach NNLL accuracy Ul9 |. 



3.3 Subleading logarithms and resummation prescriptions 

In the previous subsection we derived all the NNLO threshold corrections for qq — > ti through 
NNLL accuracy. It is interesting to study the effect of subleading [1/54]+ and S(s4) terms that 
come about when inverting the cross section from moment to momentum space. As we will see, 
this is intimately related to the disagreements between various resummation prescriptions that 
have been proposed. These prescriptions are needed to avoid the soft gluon divergences in the 
resummed cross section that appear when a s reaches the Landau pole. As we have discussed 
before, in a finite-order expansion, as presented in this paper, there are no divergences and the 
results are prescription independent. 

There are three resummation prescriptions available in the literature. The earliest is the 
x-space formalism of Refs. |3], |3l| . The resummation was performed at leading logarithmic (LL) 
accuracy in momentum space and a cutoff was chosen to avoid the divergence. In practice, the 
cutoff was chosen so that numerically the resummed result would agree with the expansion of 
the resummed cross section through NNLO. This approach was also used at NLL accuracy in 
Refs. 0, H |12[ . The finite-order expansion is essentially the same whether the resummation 
is performed in momentum or moment space so, although the approach is quite different, in 
practice the numerical results from this approach are not inconsistent with the ones we are 
presenting here at LL and NLL accuracy. 

The second prescription is principal value resummation, originally developed for Drell-Yan 
production in Ref. f32|. A principal value prescription is used to bypass the Landau pole. This 
approach was used at LL accuracy for top quark production in Refs. [|, [L8|. Numerically the 
results are similar to those of Ref. . This approach has not yet been used at NLL accuracy 
for top quark production at present. 

The third approach is the minimal prescription of Ref. ||. It has been applied at both 
LL || and NLL accuracy |10| to heavy quark production. Numerically it differs substantially 
from the results of Refs. [||, §]. As discussed in Refs. |5|, [L8| this difference emerges from extra 
subleading terms, which are kept in the minimal prescription approach, that come from the 
inversion of the resummed cross section from moment to momentum space. 

Let us begin our study of subleading terms by rewriting the MS NLO corrections for qq — > 
QQ in Eq. (|3.1|) in the shorthand notation 



^ (1) (s 4 ) 
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«.s- 



7T 



ci<5(s 4 ) + c 2 
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s 4 



+ c 3 



ln(s4/m 2 
s 4 



(3.14) 



with c 3 = 4C F , c 2 = 2ReT'H - 2C F + 2C F \n(sm 2 / (t lUl )) - 2C F ln(/4/m 2 ), and c x 
This result actually comes from the inversion of the moment space expression 



^(l)ggS+V 
'MS 



£(D(iV) = a B^ + cMN) + C3 J l(iV ) } = a B 
71 



TT 



Cl -c 2 lnA>+y(ln 2 iV + C2; 



, (3.15) 



where iV = iVe 7J5 , with 7^ the Euler constant, and Jo and I± are given in Appendix A. Eq. 
( p. 15] ) comes directly from the expansion of the resummed cross section in moment space in 
Eq. 
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Now, let us examine the NLO expansion at NLL accuracy, with /j,f — \ir — m. At that 
accuracy, we keep only the In 2 N and In N terms in Eq. Q3.15 ). Since IniV = IniV — je, upon 



inversion to momentum space we get back the [In(s4/m 2 )/s4] + and [1/54]+ terms in Eq. ( ft.l4| ) 
plus the following extra terms: 



a — 

71 



C37l , C2 



+ C27B - —C 3 



5(s 4 ) . (3.16) 



But there are no terms involving 7^ in the exact NLO calculation, i.e. in the term c\ in Eq. 
( |3.14| ) which comes from cr^ 95S+V ; therefore these terms are clearly unphysical. They should 
not appear in the cross section because of the definition of the MS scheme. The je terms are 
an artifact of the inversion from moment to momentum space. Indeed, if we had kept NLL 
terms in IniV rather than IniV, there would be no je terms. Also the coefficient of the (2 
term is wrong. As can be seen from the full NLO corrections, it has the wrong sign. We can 
study the numerical effect of these extra terms on the cross section. At the Tevatron, with 
y/S = 1.8 TeV and m = 175 GeV/c 2 , the full NLO S + V corrections, Eq. (Q, are 1.26 pb 
for Hf = = Tn, see Table 1. At NLL accuracy, the corrections are 1.31 pb. If we include 
all the unphysical terms of Eq. fl3.16j ), the NLL corrections become 0.39 pb, clearly very far 
from the true size of the full S + V corrections. Keeping only the (2 term in Eq. ( |3.16| ) we 
find 0.74 pb, which is closer to but still well below the full S + V corrections. If we keep 
only the (2 term but with the opposite sign, as in the full corrections, we find 1.88 pb. It is 
clear that keeping unphysical subleading terms in the cross section can produce very misleading 
results. Even if we discard the unphysical terms and keep only some of the physical subleading 
terms we can still make erroneous predictions, especially if the coefficients are wrong. One 
of the greatest achievements (or, from another viewpoint, tests) of the formalism of theshold 
resummation at NLL and higher accuracy is that it accurately reproduces the exact NLO cross 
section both analytically and numerically. In fact, one may argue that only because of this 
agreement is threshold resummation worthwhile. After all, if the corrections to be resummed 
are not dominant then the necessity for resummation is greatly diminished. That holds not 
only for heavy quark production but also for other QCD processes such as direct photon JZIJ, 
W+jet pPf , and single-jet production |p4}] . Therefore it is important to ensure that we don't 



introduce terms, unphysical or otherwise, that would spoil this agreement. 

We can extend our study of subleading logarithms to NNLO. The NNLO threshold correc- 
tions in moment space are given in shorthand notation by 



a^(N) = a B ^-\^N-c 3 c 2 ln 6 N + 
2ir z 4 



ln 2 iV 



2c 2 lniY(c 1 + ic2C 3 ) + U + ^Cs) + T% K, 2 - loop)! ,(3.17) 



i.e. by the square of the terms in curly brackets in Eq. (|3.15|) , plus a function F that com- 
prises the (3q terms that come from changing the argument in the running coupling, a s (/i' 2 ) = 
a s (n 2 )[l — /3 ln(/i' 2 //i 2 )(a; s (/^ 2 )/(47r)]; the two-loop K terms, with K defined below Eq. fl2~8[ ); 
square terms from the off-diagonal soft anomalous dimension matrix elements; and two-loop 
Tg and other terms (we have also absorbed in F a term — T 2 , with T\ = ci(/i = m)). Note 
that apart from the K terms, the other two-loop terms appear only beyond NNLL accuracy 
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and are not known at present. We then rewrite Eq. ( 3.17|) in terms of ^3, 12, h, and Iq denned 
in Appendix A. 

We can then immediately invert back to momentum space and find 
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(3.18) 

Here F denotes the terms that come from the inversion of F and starts contributing at NLL and 
higher accuracy. One can easily see that this form, with the appropriate explicit expression for 
F, agrees with the NNLO-NNLL expansion given in Eq. (|3.8|) , including the additional scale 
terms of Eq. (|3.11|) . Clearly we don't know all the [l/s 4 ] + and 5(s 4 ) terms because of unknown 
two-loop corrections in F. Of course there are no 7^ terms in the NNLO cross section. This 
is also known from the two-loop Drell-Yan cross section |33], B3L BSfl: many of the NNLO terms 



are the same as for top quark production since the exponent in the resummed cross section that 
comes from the incoming partons is universal; and there are no 7^ terms. Again, this follows 
from the definition of the MS scheme. These terms are clearly unphysical at both NLO and 
NNLO and indeed at any higher order. 

Now, let us see what happens if one keeps the logarithms only at a certain accuracy. At 
NNLL accuracy, we keep the ln 4 iV, ln 3 iV, and In 2 N terms in a (2) (iV), Eq. ( [OTP . Then, the 
subleading terms from the inversion to momentum space are: 
,2 
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The subleading ( terms in Eq. Q3.19 ) appear also in the full cross section of Eq. (|3.18| ). 
However, a comparison between these two equations shows that some of these terms have the 
wrong coefficients. The terms in Eq. ( |3.21|) are again the unphysical 7^ terms that should not 
appear at any order of the perturbative series. Again, if we had kept NNLL terms in In N rather 
than In AT, there would be no je terms. The (3q and K terms in Eq. ( |3.20D would also be absent 
in an exact calculation. This is because they appear in integrals of the form of Eq. (A.l), which 
upon inversion to momentum space should give back the original "plus" distributions with no 
subleading terms. 

We can again study the numerical effect of these extra subleading terms. In the following, 
we keep fip = fiR = m. At the Tevatron, with v^S* = 1.8 TeV and m = 175 GeV/c 2 , the 
MS NNLO-NNLL corrections for qq — > ti are 0.80 pb as we saw in Table 1. If we keep also 
the subleading ( terms in Eq. ( |3.19| ) the result becomes 0.39 pb. If we keep the subleading 
terms in both Eqs. ( |3.19|) and ( |3.20|) , the corrections become 0.13 pb. Finally, if we include 
all subleading terms, Eqs. ( |3.19| ), ( |3.20| ), and ( |3.21| ) the corrections become 0.08 pb. This last 



result is similar to the result presented in Ref. [10| (note that different parton densities are used; 
also our formalism resums the fully differential cross section while [TO] resums only the total 
cross section; the latter approach introduces some additional errors, see the discussion in Ref. 
|36|| ). Clearly the inclusion of the unphysical 7^ and other terms decreases the NNLO-NNLL 
corrections by a factor of ten. The effects of these unphysical terms are much bigger than 
those of the LL, NLL, and NNLL terms. It is difficult to accept a result in which unphysical 
subleading terms dominate the three leading powers of the logarithms. Such a result defies the 
meaning of leading level, next-to-leading level and so on. And as was evidenced by the NLO 
exercise, a result with these subleading terms substantially underestimates the correct value 
for the cross section. We also note that if we keep the subleading ( terms from the inversion 
at full accuracy, as in Eq. ( |3.18j ), the corrections are 0.45 pb, much closer to the NNLO-NNLL 
result. Of course, we can't derive the full NNLO cross section beyond NNLL accuracy because 
of missing two-loop terms, but this certainly indicates that the corrections with subleading 
terms tend to get larger the better the accuracy. 

We can also repeat this exercise at NLL accuracy. Here we will disregard the 7^ and (3 , K 
terms. At NLL accuracy for a^ 2 \ the subleading ( terms from the inversion are 



-C3C2C2 



(3.22) 

Again, the subleading terms above appear also in the full cross section of Eq. ( |3.18| ) but some of 
these terms have the wrong coefficients. As we saw in Table 1, the NNLO-NLL corrections are 
1.01 pb. If one adds the subleading terms above, the result becomes 0.23 pb. Our conclusions 
remain the same. We also note that if we keep only the [(ln(s 4 /m 2 ))/s 4 ] + subleading terms 
from the inversion, we find that the corrections are 0.58 pb, quite different from the result we 
get (0.80 pb) when we calculate the full NNLL terms (if we also include the unphysical terms 
the disagreement is far worse). In addition, we observe that if we perform a LL calculation and 
then include subleading NLL (i.e. [(In 2 (s4/m 2 ))/s4] + ) terms from the inversion involving je, 
the corrections become 0.28 pb versus the 1.01 pb that we find with a full NLL calculation. All 
these exercises highlight the numerical problems that are encountered if one includes unjustified 
subleading terms in the expansion. 
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As is noted in Ref. [37], the integrals of the leading logarithmic distributions with any 



smooth function, such as the convolution with parton distributions, produce factorial contribu- 
tions at nth order of the form a™ (2n — l)\/n\ + .... These factorial terms naturally arise in both 
the exact cross section at any fixed order and in the finite-order expansions of the resummed 
cross section. In the minimal prescription of Ref. || subleading terms are kept in the resummed 
cross section and its expansion in order to avoid certain power corrections, arising from these 
factorial contributions, which have been shown to be absent in the Drell Yan cross section in 
Refs. 0, |39fl. This is indeed a problem that deserves attention as pointed out in Ref. ||. 



However, the absence of unphysical power corrections does not require the introduction in the 
expansion of unphysical terms, which, as we have seen above, may greatly underestimate the 
true value of the cross section. The exact cross section does not have these power corrections 
but also it does not have the je terms, only the ( terms, as we have seen explicitly at NLO 
and NNLO. Moreover, one has to be careful not to introduce extra terms, even physical ones, 
that produce erroneous numerical results. At low orders, with n — 1,2 the factorial contribu- 
tions are negligible or small anyway, and certainly smaller than other terms in the expansion 
(see also discussion in Ref. []1&|] ). At higher orders of course the factorial contribution grows 
and moreover we have ever increasing numbers of uknown coefficients of subleading logarithms 
(which actually may be more important numerically than the factorial terms). Therefore we 
stop the expansion at NNLO, thus avoiding the theoretical problem with power corrections. 
At NLO and NNLO we can trust the perturbative expansion, as we have determined all log- 
arithmic coefficients. Our numerical analyses above confirm that. We also note that in our 
approach the results for the NNLO corrections do not change substantially when going from 
LL to NNLL accuracy (they change from 0.59 pb to 0.80 pb) while the corresponding results 
with all subleading terms change by almost an order of magnitude (that can also be seen by 
comparing Table 1 of Ref. |J with Table 2 of Ref. flit)] ). The relative stability of our results 
versus logarithmic accuracy is an additional justification of our approach. We can investigate 
keeping the ( terms as a rough estimate of error, as we have done above, but we should keep in 
mind that they may not necessarily offer an improvement on the calculation as evidenced from 
the NLO results. Therefore, we find it best in the numerical analyses presented in this paper 
not to retain any terms beyond NNLL accuracy. Thus, we do not find the very fast convergence 
of the higher order corrections that is claimed in Refs. || |l(| . 

In a recent paper [[HJ it is argued that threshold enhancements are dominated by the region 
where the hierarchy among different powers of the threshold logarithms is lost, and therefore 
NLL resummation is not reliable. Our numerical results at fixed order do not agree with this 
claim, although it is certainly true that the coefficients of lower powers of the logarithm can be 
large. Moreover, we note that even if the hierarchy among different powers of the logarithms 
were lost, at NNLO we have determined the coefficients of all the powers of the logarithms, 
so our results are reliable regardless. Beyond NNLO, however, there are subleading powers of 
logarithms with undetermined coefficients, which can be large, and then the ambiguities with 
regard to the effect of subleading terms increase. Therefore, as discussed above, for detailed 
numerical results we prefer to stop the expansion at NNLO. 
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3.4 NNLO-NNLL threshold corrections for gg — > QQ 

For the gg channel in the MS scheme the NNLO-NNLL corrections are 
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(3.23) 



where C|^f 9S+V * s defined in analogy to Eq. ( |3.9|) , and the elements of the soft anomalous 
dimension matrix (dropping gauge- dependent terms and an overall a s /7r) are 
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As for the qq channel, we can also derive at NNLL accuracy additional [l/sj 
involving the scale. These terms are 



and 5(54) terms 
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where T-^ 99 is defined in analogy to its counterpart in the gg channel. 

At the Tevatron, with \/£ = 1.8 TeV and /ip = [1r = m = 175 GeV/c 2 , the exact NLO cross 
section for the gg channel is 0.55 pb and the NNLO-NNLL corrections provide an additional 
0.32 pb. The relative size of the NNLO corrections for the gg channel compared to NLO is 
much greater than for the qq channel. This is because of the different color coefficients in the 
expressions for the two channels, as is obvious from the coefficients of the leading logarithms. 

NNLO results for the gg channel in PIM kinematics are presented in Appendix C (see also 
Ref. ill). 



3.5 Top production at the Tevatron 

In this subsection we add the numerical contributions from the qq and gg partonic channels and 
present some numerical results for the top quark total cross section and transverse momentum 
distributions at the Tevatron (see Appendix B for a discussion of the hadronic calculation). We 
use the CTEQ5M parton densities |3(| . 



In Fig. 3 we plot the Born term and the NLO and NNLO corrections for top quark pro- 
duction at the Tevatron with = 1.8 TeV as a function of the top mass. At NLO we show 
the exact corrections as well as the NLL threshold corrections and the full S+V threshold 
corrections. At NNLO we show results with both NLL and NNLL accuracy. 

In Fig. 4 we plot the exact NLO and the NNLO-NNLL MS top quark cross section at the 
Tevatron with \/S = 1.8 TeV as a function of the top mass. We note a dramatic decrease of 
the scale dependence of the cross section when we include the NNLO-NNLL corrections. We 
observe that the NNLO cross section is uniformly above the NLO cross section for all scale 
choices (we stress that there is no field-theoretical reason for the NNLO results to lie within the 
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Fig. 3. The Born, NLO, and NNLO corrections for top quark production at the Tevatron with 
V5= 1.8 TeV. 
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Fig. 4. The total cross section for top quark production at the Tevatron with =1.8 TeV. 
The labels are as in Fig. 1. Recent results from CDF and DO are also shown. 
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Fig. 5. The total cross section for top quark production at the Tevatron with \/S = 2.0 TeV. 
The labels are as in Fig. 1. 



NLO results). We also show recent results from CDF JIT] and DO and note the agreement 
between experiment and theory. In Fig. 5 we show the corresponding results for the upgraded 
Tevatron with v^S" = 2.0 TeV. 

In Table 2 we list the exact NLO and the NNLO-NNLL total cross sections in pb for top 
quark production at the Tevatron with y/S = 1.8 TeV and 2.0 TeV, a top mass m = 175 
GeV/c 2 , and scale [i = m, m/2, 2m. The NNLO-NNLL cross section with \fS = 1.8 TeV is 
6.3 pb versus 5.2 pb at NLO, an enhancement of over 20%, at // = m. Good agreement is 

in 



176.1 ±6.6 GeV/c 2 



and 



observed with recent results from CDF: a = 6.5_!; 4 pb with m 
DO: a = 5.9 ± 1.7 pb with m = 172.1 ± 7.1 GeV/c 2 |[42|| . Similar enhancements are noted for 
the upgraded Tevatron energy. 



pp — > tt 


v / 5= 1.8 TeV 


= 2.0 TeV 




NLO NNLO 


NLO NNLO 


fi = m/2 


5.4 6.4 


7.4 8.9 


fi = m 


5.2 6.3 


7.1 8.8 


fi = 2m 


4.7 6.3 


6.5 8.8 



Table 2: The hadronic tt production cross section in pb for pp collisions with \/S 
and 2.0 TeV, and m = 175 GeV/c 2 . 



1.8 TeV 



We would like to stress that the significantly reduced scale dependence should not be in- 
terpreted as an equivalent reduction of the uncertainty in the value of the cross section. Other 
sources of error such as from parton distributions, subleading logarithms, and distance from 
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Fig. 6. Top quark transverse momentum distribution at the Tevatron, with \/S = 1.8 TeV and 
2.0 TeV, and m = 175 GeV/c 2 . The labels are as in Fig. 1. 



threshold, can provide more uncertainty than the scale variation; and those errors cannot all be 
calculated precisely at present. This is why we do not give more than one significant figure after 
the decimal point in the numbers of Table 2. However, it is gratifying to see that perturbation 
theory behaves as we would expect it to H3[: at higher orders the scale variation decreases. 



Since the effect of subleading logarithms is the greatest uncertainty in the calculation, the total 
ti cross section can be written with a minimal error estimate as 6.3ig4 pb at v^S* = 1.8 TeV 
and as 8.8ta5 pb at v^S* = 2.0 TeV, where the larger lower error indicates the uncertainty from 
subleading terms. 

Our formalism allows the calculation of any relevant differential cross section. Transverse 
momentum and rapidity distributions with leading logarithmic resummation have been pre- 
sented for top production at the Tevatron in Ref. ||44|| . The exact NLO and the NNLO-NNLL 
top quark transverse momentum {p T = t\U\j 1 s 2 — m 2 ) distributions at the Tevatron, with 
\fS =1.8 TeV and 2.0 TeV, and m = 175 GeV/c 2 are shown in Fig. 6, again in the MS 
scheme. We note an overall enhancement at NNLO with little change of shape. Similar con- 
clusions are also reached for the rapidity distributions E5 . 



Finally, we note that threshold resummation is also relevant for bottom quark production 
at the HERA-B experiment. Leading logarithmic resummed results for the bottom quark 



total cross section and differential distributions have been presented in Refs. p7| , [46[j ; for the 
NLL resummed cross section see Refs. [|, 0. At NNLO-NNLL with y/S = 41.6 GeV and 
fi = m = 4.75 GeV/c 2 we find a cross section for that experiment of 42 nb, while the NLO 
cross section is 18 nb. 



4 N 3 LO threshold corrections 

We now go beyond the NNLO corrections and expand the resummed cross section to next-to- 
next-to-next-to-leading order at NNLL accuracy. 
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For the qq channel in the MS scheme the N 3 LO-NNLL threshold corrections are given by 
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(4.1) 



where ci, c 2 are defined below Eq. ( |3.14j ). Note that at NNLL accuracy there are no cubic terms 
in Ts; they start contributing at 0([(ln 2 (s 4 )/m 2 )/s 4 ] + ). One can also derive terms involving 
the factorization and renormalization scales at lower powers of the logarithms as was explained 
in Section 3. 

For the gg channel the leading logarithms at N 3 LO are of course the same as for the qq 
channel with the substitution C F —>■ Ca in the coefficients. Beyond leading logarithms the 
more complex color structure of the hard scattering for the gg channel makes the calculation 
more lengthy, as is evident already at NNLO, and will not be pursued here. 

In the DIS scheme the corresponding result is 
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(4.2) 



where c[ = c^}f S+Y , and c' 2 = 2ReT'g - C F /2 + C F \n(s 2 / (hm)) - 2C F ln(/4/m 2 ). 

We can extend our study of subleading logarithms to N 3 LO. The N 3 LO threshold MS 
corrections for qq — > QQ in moment space are given in shorthand notation by 



a 
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6tt 3 



ci -c 2 lniV + ^(ln 2 iV + C 2 ) 
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(4.3) 



i.e. by the cube of the terms in curly brackets in Eq. ( |3.15| ) plus a function F' that gives 
the /3q terms that come from changing the argument in the running coupling, the two-loop K 
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terms, additional square and cubic terms in the soft anomalous dimension matrix elements, 
and two- and three-loop T$ and other terms. Note that we have also absorbed in F' the terms 
2T 3 - 32? [c! - c 2 In JV + (c 3 /2)(ln 2 N + C2)], with T x =c 1 (fi = m). We then rewrite Eq. flC 
in terms of J5, 14, 13, 12, Ji, and Iq defined in Appendix A. 
Inverting back to momentum space, we have 
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where F' comes from the inversion of F' and starts contributing at NLL and higher accuracy. 
One can easily see that the above equation, with the appropriate explicit expression for F', 
agrees with the N 3 LO-NNLL expansion given in Eq. (44.). Of course we cannot derive all of 
the [(ln 2 (s 4 /m 2 ))/s 4 ] + and lower terms in Eq. (|4.4p because of unknown 3-loop corrections in 
F'. 

Again, let us see what happens if one keeps the logarithms only at a certain accuracy. At 
NLL accuracy for the subleading terms from the inversion (keeping only £ terms and no 
Po,K,je terms as discussed in Section 3.3) are: 
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2 6 \ 8 3 4 

At NNLL accuracy for the corresponding subleading terms from the inversion to mo- 
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mentum space are: 
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A comparison of both the NLL and the NNLL subleading terms with Eq. (|4.4| ) shows that 
most of these terms have incorrect coefficients. The subleading terms bring down the value of 
the MS N 3 LO corrections for qq — > tt at the Tevatron with \J~S = 1.8 TeV and m = 175 GeV/c 2 
from 0.9 pb (1.1 pb) at NNLL (NLL) to around 0.3 pb. As we discussed in Section 3.3, retaining 
subleading terms with incorrect coefficients in the expansions can produce misleading results. 
We also note that in the DIS scheme the corresponding corrections (without subleading terms) 
are 0.2 pb, again smaller than the MS result because of the specification of the DIS scheme. 

Since we don't know the (potentially large) coefficients of subleading powers (beyond NNLL) 
of the logarithms at N 3 LO, in contrast to the NNLO calculation where all logarithms were 
determined, we cannot make firm numerical predictions at this order. This also relates to 
the questions raised in Ref. 0] as we discussed earlier. A full N 3 LO calculation may give 
substantially different (and smaller) results from the N 3 LO-NNLL calculation; but at present 
we cannot calculate corrections beyond NNLL accuracy. Therefore, as we have stated before, 
for reliable numerical predictions we prefer to stop the expansion at NNLO. 



5 N 4 LO and higher-order threshold corrections 

Finally, we briefly discuss the corrections at next-to-next-to-next-to-next-to-leading and higher 
orders. 

For the qq channel in the MS scheme at NLL accuracy, the N 4 LO threshold corrections are 
given by 
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By matching with the exact NLO cross section, as we have described before, one can derive the 
full NNLL terms as well. We note that no cubic or quartic terms in Ts appear at NNLL accuracy. 
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A full determination of the cross section at this order would require four-loop calculations. The 
leading logarithms for the gg channel at N 4 LO again follow from Eq. (|5.1|) with the substitution 
Cf — * Ca- 

In the DIS scheme the corresponding NLL result is 
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Again, numerically the corrections in the DIS scheme are smaller than in the MS scheme. 

The finite-order expansion procedure can be extended to arbitrarily high orders at NNLL 
accuracy. In general, at nth order in a s (beyond the Born term) the leading logarithms in the 
MS scheme are 
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For the DIS scheme we only need to replace (2Cp) n by C F in the above equation. For the gg 
channel we simply replace Cp by Ca as discussed before. It is easy to check that Eq. (|5.3| ) 
reproduces the leading logarithms in all the expansions presented in this paper. 



6 Conclusions 

Threshold resummation can make powerful improvements to NLO calculations of heavy quark 
production cross sections. The analytical form of the threshold corrections to the fully differen- 
tial cross section can be derived at NNLL accuracy at arbitrarily high order and explicit results 
have been provided in this paper through N 4 LO. The role of subleading logarithms has been 
studied in detail and it has been shown that care must be taken to arrive at reliable numeri- 
cal predictions for the cross section. For top quark production at the Tevatron NNLO-NNLL 
predictions have been made for the total cross section and transverse momentum distributions. 
The NNLO-NNLL corrections are significant and they dramatically reduce the dependence of 
the cross section on the scale relative to NLO. 
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A Mellin transforms 

Here we present some useful formulas for the Mellin transforms that are used in the resummed 
cross section and the finite-order expansions. 
We define 
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The results below are identical for either definition. 
Then, we have [Ef 
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Here N = Ne" fE with >y E the Euler constant, 7^=0.577... Also C2 = vr 2 /6, (4 = vr 4 /90, Ce 
tt 6 /945, Cs = vr 8 /9450, etc., while Cs = 1.2020569..., Cs = 1.0369278..., (7 = 1.0083493..., etc. 
Then 
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The expressions for I n have been presented up to n = 3 in f2|. Here we extend this table to 
n = 7, useful through N 4 LO expansions: 
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B Hadronic cross sections 



The double differential hadronic cross section d 2 Oh x } l2 j dT\d\J\ is written as a convolution of 
parton distributions with the partonic differential cross section: 



S 



2 (Puh^iSi Ti, U\) 



dTi dUi 



E 



'■j 



1 dXl } dx 2 



X\ 



x 2 



. 2 \ i 1 2 \ 2 d &ij( S > ^1, Ui) 



dt\ du\ 



(B.l) 

where the sum is over all massless parton flavors and (p^x,^) are the parton distribution 
functions for flavor i in hadron h at factorization scale p,p. The hadronic invariants S, T 1; U\ are 
defined in analogy to their partonic counterparts. The lower limits of the momentum fractions 
of the partons in the hadrons are given by 27 = —U\/(S + Ti) and x 2 = — X\T\f(x\S + U\). 

By making a transformation from the variables (Ti, Ui, Xi,x 2 ) to the variables (Y, p\, x±, s 4 ), 
with Y and px the rapidity and transverse momentum, via 
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Now, let us write the /cth-order partonic threshold corrections in the shorthand notation 
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By substituting the above expression for the partonic threshold corrections in Eq. ( |B.3| ) 
we can write the hadronic fcth-order corrections as 
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After some rearrangements of terms and partial integrations, we can rewrite this as 
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The transverse momentum distributions are then given by 
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4(Pt + m2 ) / The total cross section can then be retrieved by integrating Eq. 



( B.7|) over p T with lower limit and upper limit J S/A — m 2 . 



C NLO and NNLO threshold corrections in PIM kine- 
matics 

In this Appendix we present results for the NLO and NNLO expansions of the resummed cross 
section in heavy-quark-pair inclusive kinematics (see also Refs. 0, 0, [19]). Here the distance 
from threshold is measured in terms of the variable z = Q 2 /s, with Q 2 the invariant mass 
squared of the heavy quark-antiquark pair, and the corresponding "plus" distributions are of 
the form [(ln fc (l - z))/(l - z)} + . 



C.l NLO threshold corrections 

In the MS scheme the NLO-NLL corrections for the qq — > QQ channel are 
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We note that in the 5(1 — z) contribution the expansion reproduces only the scale-dependent 
terms. The rest of the 5(1 — z) terms can only be obtained by matching to the exact NLO cross 
section in PIM kinematics ]19[] ; this is exactly analogous to what was presented in Section 3 for 
1PI kinematics. 
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In the DIS scheme the corresponding result is 
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For the gg — > channel the NLO-NLL corrections in the MS scheme are 
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C.2 NNLO threshold corrections 

We now present the NNLO corrections in PIM kinematics. We give explicit results at NLL 
accuracy. As noted in Section 3, to reach NNLL accuracy we need to derive matching terms 



in PIM kinematics [19]. The only other difference between the expansions in the two different 
kinematics are in the extra terms involving ]n(t\U\/m A ) for the 1PI fixed-order expansions 
relative to the PIM expansions. 

In the MS scheme, the NNLO-NLL corrections for qq — > QQ are 
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We note that, at NLL accuracy, in the [ln(l — z)/{l — z)) + coefficient we derive only the terms 
involving the factorization and renormalization scales, and in the [1/(1 — z)] + coefficient we can 
only determine quadratic terms in the scale logarithms. 
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In the DIS scheme the corresponding corrections are 
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In the MS scheme for the gg — > QQ channel the NNLO-NLL corrections are 
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